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1.0 


Introduction 


This report presents the results of a study on 
video data compression techniques applicable to space 
flight communication. This study is exclusively directed 
towards monochrome (i.e. black and white) picture communi- 
cation with special emphasis on feasibility of hardware 
implementation. The primary factors for such a communi- 
cation system, in space flight applications are: 

1) picture quality, 

2j) system reliability, 

3) power consumption, 

4) hardware weight. 

In terms of hardware implementation, these are directly 
related to hardware complexity, effectiveness of the hard- 
ware algorithm, immunity of the source code to channel noise, 
and data transmission rate (or transmission bandwidth) , 

This report will recommend a system and summarize its hard- 
ware requirement. In addition, this report will provide 
sufficient data on various parameters involved. 

Simulations of the study were performed on the im- 

. . 

proved LIM Video Controller. The LIM Video Controller is 

- . 

computer-controlled by the META-4 CPU. The functional block 
diagram of the LIM Video Controller is illustrated in Figure 
1. The LIM Video Controller processes video signals recorded 
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on the Ampex Video Recorder, Model DR-10. This is done 
by first A-to-D converting the video signals (in groups 
of four horizontal lines) , transporting the digital data 
to the META-4 main storage for processing, then the pro- 
cessed signal is reconstructed by D-to-A conversion and 
finally recorded on the video disc recorder for visual 
display. The DR- 10 has 600 tracks (one frame per tiack) . 

Half of the disc allocation (i.e. 300 frames of video 
signals) is used for storage of reference signals. The 
remaining half is used for recording the processed signals. 
This corresponds to 10 seconds of real-time video signals. 

This seems sufficient for demonstration purposes for the 
effectiveness of the algorithm under study. By preserving 
the pre-recorded reference signals on the video disc re- 
corder, effectiveness of various schemes can be compared 
fairly by recording the processed signals on a video tape 
recorder. The improved LIM Video Controller and its new 
supporting software enable 300 frames of video signals to be 
processed in approximately 4-1/2 hours. Each processed frame 
has a superimposed title and reference frame number for 
ease of identification. Noise degradation due to simulation 
is minimal. However, the video tape recorder (Sony model) 
has limited bandwidth and causes certain noise effects. For a 
detailed description of the LIM Video Controller, see the LIM 
Video Controller Operations Manual.. 
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Figure 2. Typical Picture Communication Block Diagram 










Most video data compression techniques are based 
upon the above facts. One of the methods, originally in- 
vestigated by Landau and Slepian [2] , uses spatial statis- 
tical correlation exclusively. The picture is first 
partitioned into small regions of subpictures. Statis- 
tical data reduction is performed on each subpicture by 
considering each subpicture as an independent random vec- 
tor. Coordinate transformation is applied to the random 
vector. The objectives of the coordinate transformation 
is to diagonalize the original covariance matrix of the 
subpicture and to produce an orthonormal basis similar 
tb the Karhunen-Loeve procedure. Bit rate reduction is 

obtained by discarding or quantizing with fewer information 

• | 

bits those components that have lower statistical 
variances. The transformation used is the Hadamard trans- 
formation, where the basis vectors corresponds to the 
row vectors of a Hadamard matrix. This method has the 
advantage of simple and fast hardware implementation, short 
delay between the real-time and the processed pictures , and 
the coding errors due to channel noise are confined to sub 

pictures (i.e. high coding reliability) . Compression ratio 

; i 

of 4 to 1 (assuming the original video signals are linearly 
digitized by an 8-bit A-to-D converter) is achievable with- 
out substantial deterioration in picture quality. However, 
this method does not use the time statistical correlation 
between adjacent frames of the video sequence. 
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Another relatively simple method is the frame-to- 
frame differencing, variable-length coding technique [3] . 

This method only utilizes the time statistical correlation 
in recognizable video sequences. Here, a reference frame 
of full information is transmitted. For the subsequent 
frames, the video signals are compared with their corres- 
ponding video signals of the preceeding frame, and only those 
differences that exceeded certain pre-chosen threshold are 
transmitted. At the receiving side, the reference is first 
reconstructed, the subsequent frames are updated by the in- 
formation received. This method basically uses the fact that 
iflost recognizable sequences have a large proportion of sta- 
tionary objects; thus, the amount of information changed 
• | 

irom frame to frame is rather small. Source coding using 
this method requires that a substantial portion of the code 
be allocated for position markers to indicate where in the 
picture the changes take place. Moreover, a rate buffer is 
required to achieve fixed rate transmission. Furthermore, 
tihis method is very sensitive to channel noise; a decoding 
error in one frame causes errors in subsequent frames. Thus, 
to improve picture reliability, reference frames must be 
transmitted every so often. 

In a previous LINKABIT video study report [4] , LINKABIT 
provided a buffer-free technique that utilizes both the spa- 
tial and time statistical correlations of recognizable scenes. 
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This method is an essential combination of the two methods 
described. Like the Landau and Slepian method, video sig- 
nals are first transformed into Hadamard coordinates. Bit 
rate reduction is first achieved utilizing the spatial sta- 
tistical correlation; a reference frame is transmitted in 
this manner. For the subsequent 3 frames, the Hadamard com- 
ponents of the new picture are compared with the correspond- 
ing components of the reference frame and only the differences 
of! a few selected components (those that have highest statis- 
tical variances) are quantized and transmitted. The process 
is repeated every four frames. This method seems capable 
of: reproducing a scene with a compression ratio of @ to 1 
while retaining recognizability of the scene. It is the 
sole objective of this report to explore improvements in this 
technique. Due to its relatively simple hardware implementa- 
tion and partial frame storage (instead of a full frame stor- 
age required in method described above) , this method seems to 
be most promising for space flight application. (Although other 
image coding methods are available, such a 2-dimension Fourier 
transformation and other related techniques, their computational 
complexities and bulk storage required for past frame infor- 
mation limit their usage to ground instruments.) 
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3.0 


Results of Study of Two-Dimensional Video Data 
Compression 

Since the technique suggested in this report uses 
the Landau and Slepian method for spatial statistical data 

i 

reduction, it is essential to examine this method and ex- 
plore possible improvements. We shall begin with a review 
of the Landau and Slepian two-dimensional transformation 
technique. 

We assume the video signal source being processed is 
regular commercial NTSC TV. Each frame of the video signal 
consists of two fields interlaced with each other. The 
video signals are digitized by an A-to-D, converter. The 
A-to-D converter must sample at a frequency above the Nyquist 
rate of the desired bandwidth, and a sufficient number of 
lj>its of information per sample is required to ensure smooth 
Video reproduction. For monochrome TV signals, sampling at 

512 samples per horizontal line and 8-bit grey level resolu- 

i 

tion seem sufficient to reproduce reasonable quality picture 
without false image contouring. 
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Each frame of video signal consists of 525 horizontal 
lines. Among these, 45 horizontal lines are used for genera- 
tion of vertical sync pulses and are blanked. Thus, each 
frame contains, at most, 480 lines of visible information. 

When digitized accordingly, we can view each frame of digi- 
tized video signals as a lattice of 480 x 512 sample points, 
x. (i = 1, 2, ..., 480; j = 1, 2, ..., 512). Each x. . 

X J X J 

has integer representation value between 0 and 255. In 
other, words, each frame of digitized video signal can be 
represented by a vector in an Euclidean space of dimension 
480 x 512. This representation is illustrated in Figure 3. 

Since spatial statistical correlations between sampling 
points are effective only for neighboring points, it is 
desirable to partition the picture into subpictures where 
spatial statistical correlations within the subpicture are 
highest. One method of achieving this is to partition the 
picture into rectangles of size m x n (m vertically and n 
horizontally) , where m is a divisor of 480 and n is a divi- 
sor of 512. For commercial TV signals, since fields 
overlap each other, the subpictures can be formed within 
each field. In this way, subpictures of field 1 overlap with 
subpictures of field 2. An attempt to form subpictures from 
both fields may not be advantageous, since a delay of at 
least one field time is required between the real-time and 
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processed video signals. Furthermore, field 1 and field 2 
alre in reality two different pictures at two different in- 
stances; consequently, for scenes with a lot of object move- 
ments, the spatial statistical correlation between neighbor 
points with different fields may not be effective. 

Choice of m and n are noinnally determined by con- 
sidering the following factors; 

1) Higher data compression ratio can be achieved 
with larger m and n. 

2) Computational complexity (consequently, hard- 
ware implementation) and time delay between real-time and 
processed signals generally increase with m and n. 

3) For subpictures of too large m and n, the spatial 
satistical correlations between furtherest sample points 
within the subpicture diminish; thus, data compression ratio 
may not be further improved by increasing m and n. 
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Experiments have shown that subpictures of size 
4x4 seem to be most efficient. Due to the asymmetry 
resulted from the effect of field interlacing, subpictures 
of size 4x8 should be also a reasonable choice. All 
experiments in this report are based upon subpictures of 
size 4x4. • 

Each subpicture is considered as an m x n random 
vector independent of other subpictures. Symbolically, 
each subpicture of video can be represented as: 

0 < x. . j< 255 

x^j - integer ... 

(3.1) 

' Each x^j represents the digitized , '-deo signal at 
coordinate (i,j) within the subpicture. The range of this 
vector is the lattice of (256) mn integral points lying 
within the mn-dimensional cube of size 255. Due to the 
spatial statistical correlations of recognizable pictures, 
not all of the vectors in the range are comprehensible to 
th;. human eye. By considering the set of subpictures extracted 
from recognizable pictures, statistical data reduction in the 
sense of least mean square error generally can be obtained, by 
using Karhunen-Loeve procedure. This procedure requires gen- 
eration of the orthonormal basis which diagonalizes the co- 



X 11 * 

. . x . 
m 

Y = 

• 

i — 1 

<M 

X 

* * x 2n 


x ml * 

x mn 


13 


variance matrix of the ensemble of subpictures, (3.1), 
ofj recognizable pictures. However, the least mean 
square error criterion is not generally suitable as a mea- 
sure of visual fidelity. Other orthonormal bases were 

! 

sought. One of the intuitive choices, which was judged 
to be superior to the Karhunen-Loeve basis, is the Hada- 
mard basis * The Hadamard basis are vectors that are the 
rows of a Hadamard matrix. 

An n x n matrix, H, of integer entries is Hadamard 
if 

T 

H » H = nl 

i »p 

Where H is the transpose of H. I is the identity matrix 
of n-dimensional vector space. A Hadamard matrix of order 
2 , k = 2, 3, ..., can be obtained recursively as the ten- 
sor product of Hadamard matrix of order 2, H£: 

(3.2) 

k l k 2 

Therefore, a subpicture of size m - 2 and n = 2 

can be coordinate transformed using a Hadamard basis; the 
j ■ ' . .. ■ 

basis vectors are obtained from the row vectors of the Hada- 

k l+ k 2 

mard matrix of order 2 (which can be obtained by the 

! 

* k^ + k 2 tensor product of H^) . All row vectors of any Hada- 
x mard matrix have component values +1 or -1. Thus, the ortho 
normal basis can be obtained by dividing each row vector by 
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(3.3) 


the constant: 


ann 


. A 


k 1+ k 2 


The first basis vector thus formed is of the form: 
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1 , 

1 . 


1 . . . 1 


(3.4) 


+ 

m x n 


The remaining basis vectors are all of the form l//mn 
multiplied by a vector in which half of the components 
have value +1 and the remaining half have component value 
-1. Let us denote the Hadamard basis vector by 


b. . 
ID 


i = 1, 2, . . . m 
3 = 1, 2 , . . .n 


Then each subpicture of the video signal can be expressed as 

m n 

Y = V* V* C.., b.., . . . (3.5) 

i=l j=l 


where 


C. . = Y • b. . 
ID ID 


(3.6) 


is the projection of Y into basis vector b^ j . In particu 

lar, the first component, C^, has the following expression: 

m n 


C,, = 


v^rin 


£ £ 

i=l 3=1 


X i j 


(3.7) 
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This component is also known as the d.c. component of the 
subpicture, for reason that it is a constant multiple of 
the sum of individual sample video signals. The remaining 
components, C^j, have the form l/fm\ multipled by the dif- 
ference of sums of half of the x^j’s. Thus, when trans- 
formed, the Hadamard component has the following range: 


0 < C 1;L < 255 /mn 
and 


(3.8) 


-255/nin/2 < £ + 255v / mn/2 (to the nearest integer) 

for i / 1 and j ^ 1 

Given a vector, Z, in the Hadamard basis, the cor- 
responding subpicture, Y, can be obtained by inverse trans- 
formation: 


Y' = 


3-D 


X 11 • • 

' x ln 

X 21 * * 

* x 2n 

x ml . . 

. xmn 

A 


• b . . 



(3.9) 


3-D 


Data rate reduction is achieved by statistical analysis 

on the covariance matrix of the transformed Hadamard components 

{ 

on the ensemble of recognizable pictures. Statistical analysis 
of ordinary recognizable pictures, using a 4 x 4 subpicture size, 
reveals that the variance of the first component, C^, is the 
highest and is in excess of 10 to 1 in ratio to ;the next' highest 
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mentation. In this report, we shall follow 
reasoning more closely in researching for an 





logarithmically. The logarithmic response explains the 
reason why human eyes are extremely sensitive to false 
image contouring and graininess of the picture. Typical 
response of human eyes with respect to the difference of 
video signals of adjacent regions can be expressed as: 

Sgn(x)v A. log (B |x| + 1) (3,10) 


where A and B are positive real constants, and x is the 
difference of the video signals between neighboring re- 
gions. (Positive sign implies brighter signal and nega- 
tive sign implies darker signal.) Based upon this reason- 
ing, the quantization outpoints for a given number of bits 
of information should be chosen accordingly, as the in- 
verse of (3.10). 


Let the number of quantized levels be N^' + Nj + 1# 
which correspond to the integer set; 


-N^ , -N^+l, .. . , -2, -1, 0, 1, 2, . . . 

N 2 -1, n 2 . 

Then the quantization levels should be chosen as; 
-B.k 

- A-^fe - 1) for k = 0, -1, . . . , -N^ 

and 

B 2 k 

+ Agfe - 1] for k = 0 , 1, . . ., N 2 


(3.11) 
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In other words the quantization error is allowed to in- 
crease exponentially if the difference of adjacent com- 
ponents is large. The formulation (3.11) also allows 
"zero" representation. Quantization without zero repre- 
sentation has the disadvantage of introducing unwanted 
"sand paper" effect, which is the inability of the cod- 
ing to reproduce smooth images. Although some authors 
advocate to remedy this situation by introducing pseudo- 

j 

random noise such tactics seem more likely to dis- 
guise the bad by worse. Quantization with "zero" repre- 
sentation is highly recommended by this report. 

The cutpoints for the quantization (3.11), can be 
chosen as the + N 2 arithmetic means of the + N 2 ad- 
jacent pairs given by (3.11). A^, Bjf A 2 and B 2 are con- 
stants which determine the graininess and the maximum 
jcepre sen table value of the quantization. Since the amount 
of possible quantization levels available for a given num- 
ber of bits is in the form of 2^, often either 

= N 2 +l or N 2 = (ssy^etric quantization) is de- 

sirable for maximum usage of information available. Asym- 
metric quantization may increase slightly in hardware com- 

l • • 

plexity. For large k (such as 5 or above) , introducing 
asymmetry may not result in much information gain. 
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3.2 


Quantization of the d.c. Component of the 
Hadamard Transformation 

Logarithmic quantizations, though intuitively sound, 
cannot be applied directly to the d.c. component "f the 
Hadamard transformation. For, the video signals of an 
arbitrary subpicture, Y, can assume any value within the 
range of the A-to-D converter, i.e. (3.1), and since the 
human eyes are relatively insensitive to relative 
brightness; it is not advantageous to bias the 
usable video region to any extent. However, the d.c. com- 
ponent can be encoded logarithmically by using DPCM tech- 
nique. Here, the d.c. component of a subpicture is selected 
as a reference, the d.c. components of the following sub- 
pictures are coded as the differences of these signals. 

Af the receiving end, the reference is first reconstructed 
and the d.c. components of the succeeding subpictures are 
reconstructed by updating their preceeding reconstructed 
components. Since a frame contains 480/m by 512/n sub- 
pictures, the DPCM technique can be applied horizontally, 
vertically or both. "Horizontally" means a reference is 
sent at the beginning of each horizontal group of subpic- 
tures. "Vertically" means a reference is sent at the be- 
ginning of each vertical group of subpictures. Only one 
reference is sent for each frame if both horizontally and 
vertically. In addition, DPCM vertically requires suffi- 
cient data storage for one line group, this is due to the 






"Si 
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fact that the TV signal is horizontally oriented. 


Encoding the d.c. component using the logarithmic 
DPCM method has the following advantages: 

1) Smooth transition of grey levels between 
adjacent subpictures, which is characterized by small 
video amplitude variations, can be achieved with fewer 
bits than would otherwise be required. The smallest 
quantum jump using quantization (3.11) is 

B 1 

(e - 1) for - (light to dark) 

R (3.12) 

( e “ f° r + (hark to light) 

2) When the grey level transition between subpic- 
tures is high, it is approximated logarithmically by the quan- 
tizer; due to the logarithmic response of human eyes, the 
corresponding error of visual sensation is relatively low. 


j 3) The logarithmic quantizer can correct a step 

! 

function to within the error given by (3.12) in M 

■ j . • . • 

steps, where M is a logarithmic function of the amplitude 
of the step function and the graininess (3.12). This enables 
fewer bits of information to be allocated for quantization 
of the d.c. component than that would otherwise require. 

1 The above statements can also be explained as follows: 
Since the d.c. component of the Hadamard transform corresponds 
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Amplitude 



to an approximation, of the original picture by subpictures 
(i.e. constant grey level) , the spatial statistical 
correlation between subpictures remains valid (if the sub- 
pictures were small enough) . This spatial statistical 
correlation is in the form of smooth transitions and log- 
arithmic error tolerance of the ‘human eyes. As in most 
cases, this spatial statistical correlation enables data 
reduction by using DPCM. 

Our experiments showed that a 5-bit DPCM coding on 
the d.c. component can. effectively reproduce the picture 
with reasonable visual quality that is free from "false 
image contouring." In contrast, at least seven bits are 
required if linear quantization were used. 

' - ■ - . . f ■■ 

DPCM coding has the disadvantage of being sensitive 
to the channel noise . This is caused by the fact that 
the reconstruction of an element at the receiving end is 
dependent on the reconstruction of the previous elements. 

DPCM! horizontally will confine the errors to within a hori- 
zontal group of subpictures, and this method should be used. 

3. 3 Quantization of non-d.c. Hadamard Components 

All non-d.c. Hadamard components are constant multiples 
of differences of two sums of' video amplitudes within the 
subpicture. These components are normally very small and 
close to zero for recognizable pictures. In terms of visual 
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response, these components become significant only if 
the subpicture contains edge or spot information. Thus, 
their presence localized (or clustered) around the boun- 
daries between distinct objects. Since edges or spots 
within subpictures over the ensemble of recognizable 
pictures are generally uncorrelated, we can assume the 
spatial statistical correlation between corresponding 
components of adjacent subpictures is insignificant, 
t,hus , coding these components with DPCM may not be ad- 
vantageous. Due to their zero mean values and small 
variance (over the ensemble of recognizable pictures) , 
quantization of these components should be chosen about 
the value zero. Quantization without zero representation 
(or nearly zero representation) will result in a "sand- 
paper" effect. Zero or near zero quantization representa- 
tion for these components is recommended. Since, given a 

I 

nonzero number of bits of information, the possible quanti- 
zation levels are always even, zero representation always 
causes uneven quantization between positive and negative 
values. For components quantized with many bits (3 or more) 
the effect of the shift is negligible, but for quantiza- 
tion with fewer bits (2 or 1) , the bias due to the shift 
may not be desirable, in such cases near zero quantization 
representation should be used. Often it is advantageous 
to combine the quantization tables of coarsely quantized 
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components with that of more finely quantized components. 

For example: 12 bits are normally used to code 2 components 

with 16 quantization levels (15 outpoints) and 2 components 
with 4 quantization levels (3 cutpoints) each. The same 
number of bits is sufficient to code 1 component with 13 
quantization levels, 1 component with 12 quantization levels 
arid 2 components with 5 quantization levels. Sharing bits 
enables most efficient use of amount of information avail- 
able for a given allocated number of information bits. 
However, using this method will increase the arithmetic 
computation and subsequently the overall hardware complexity 
3.4. Selection of Quantization Levels 


Based upon the above discussion, quantization levels 
for experiments performed in this report are chosen accord- 
ing to formulation (3.11). First the number of quantization 

K 

levels is determined. This is (for non-sharing case) 2 , 
where. K is the number bits allocated. That is 


2 1 N i + N 2 + 1 


The minimum quantum jump, i.e., the graininess, using 
this method is given by (3.12) and the maximum change is 
given by 


A^ 

B 1 N 1 
(e 1 1 

- 1) 

for 

A 2 

, B 2 N 2 
(e 

- 1) 

for 


(3.13) 




■&£z*4WhXr^t&-. 





and = B£ are chosen for symmetry 
between hard and light directed transitions. The maxi- 
mum change and the graininess determine the constants A 
and B and vice versa. The maximum change for the d, c. 
component is determined by the range of i.e. 

0 <_ <_ 255 v / mn 


Normally, A^ = 


The maximum change is selected as a fraction of 
255>/mn. Here we used the fact that the probability for 
the transition of two adjacent subpictures from extreme 
darkness to extreme brightness is small. For the non d.c. 
components/ the maximum change of quantization is chosen as 
a fraction of the range of these components/ i.e., a frac- 
tion of 255/mn/2. The fraction is generally determined by 
the number of quantization levels allocated and judged 
solely by experiments. 

t 

Applying the above discussion to 4 x 4 Hadamard 
transform./ we come up with the following scheme. The com- 
ponent designation is shown in Figure 6. Here 


has 32 quantization levels using 5 bit DPCM. 
C 12 and C ^ have 7 quantization levels. 

t . . • . ■ • 

and C 31 have 15 quantization levels, 
and have 9 quantization levels. 


(7 levels) 


(15 levels) 


(9 levels) 



C 21 C 22 C 23 


(7 levels) (discard) (discard) 



C 41 C 42 C 43 

(9 levels) (discard) (5 levels) 



(discard) 


Figure 6. Hadamard Component Designation 
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C 12 ' C 13' C 14' (likewise C 21 , C 31 , and C 41 ) 
share 10 information bits. 

C 33 , C 34 , and C 43 each have 5 quantization 
levels, and the } 1 share 7 information bits. 

The overall bit requirement is 
5 + 10 + 10 + 7 = 32 bits 

or 2 bits per picture element. 

The quantization levels and their outpoints are 
shown in Table 1. 

3.5 Further Possible Improvement in Coding Efficiency 

In the above discussion, we have assumed that each 
Hadamard component can have independent occurrence within 
their range (3.7) and (3.8). To reproduce the true, ori- 
ginal picture, each component must be capable of covering 
the entire range. Source encoding by truncating the range 
always results in degradation of the reproduced picture, 
although the statistical data reduction method intends to 
limit the degradation to areas that occur rarely. Yet, 
such effects are generally felt by an observer. Hence, 
picture quality improvement can be achieved if the true 
ranges of each Hadamard component can be established. Since 
the subpicture, Y, has range: 

( X ll“* X ln\ 0 < x ij < 255 

X nl * * * X mn ) x ij ” inte ^ er 
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Table 1. Quantization Table 
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c ll ! DPCM Logarithmically by 5 Bits (Continued) 
Outpoints Representative Value 

-147 
-199 
-209 
-363 
-488 
-655 
-877 

Non-d.c. components have range: -512 to 

C^ 3 and C 31 : Quantized by 15 Levels 

Outpoin ts Representative Value 



+150 

+122 




^ +94 

+76 


! | 


+59 

+47 



i 

+35 

+28 




+20 

+15 




+10 

+7 




+4 

+2 





0 


Table 1. Quantization Table (Continued) 


-124 

-169 

-229 

-310 

-417 

-560 

-750 

■1000 
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which lies within the mn-dimensional cube of size 255. 

When transformed into the Hadamard components, the cor- 
responding boundary condition must also be satisfied. 

It is not difficult to see that the transformed Hadamard 
components satisfy the following boundary condition 

m n ___ 

53 i c iji~i c iii + i c ii"^^i - 255 ^ /2 • * • (3 * 14 

i=l j=l f 

See Appendix A. (3.14) merely states that if some com- 
ponents have very high absolute values (i.e. , near the 
boundary of representable video signals) , then the remain- 
ing components cannot have large absolute values. In other 
words, using (3.14) , we can improve the estimation for the 
ranges of occurrence of non-d.c. components of the Hadamard 
transformation by disregarding the nonaccessible region. 

To apply this method, first we have to order the 
priority of the Hadamard components. This can be normally 
chosen according to the variances. For 4 x 4 Hadamard 
transforms, we can choose the following order: 




C^ 3 to A^ = - 255v / ran/2[ 


C 31 as a ratio of: 


to A2 — A^ — [ Q (C 13 ) 


when Q(C^ 3 ) = inverse quantized representation of C^ 3 


as a ratio of: 


c i4 to A 3 _ A 2 ~ I Q ^ c i4 ) I 

when Q(C^) = inverse quantized representation of C^, 


and etc. 


Alternatively, given a set of quantization, a cor 
rection factor based upon A^, A.^, etc., is premultiplied 
by C^ 3 , C^, etc., and. these are compared with the pre- 
chosen cutpoints as in the usual manner. The inverse 
quantization representations are obtained by the repre- 
sentation value and the reciprocal of the correction fac 
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Buffer-Free Frame— to-Frame Data Compression 


This technique, proposed in a previous LINKABIT 
video study report [4], utilizes both the spatial and 
time statistical correlations of the emsemble of recog- 
nizable pictures. In this method, a reference of the 2- 
dimensional compressed Hadamard components is trans- 
mitted; for the subsequent k (k is a positive integer) 
frames, only the differences between a few major com- 
ponents and their corresponding components of the ref- 

V 

erence frame are quantized and transmitted. The process 
is then repeated every k + 1 frames, k was chosen to be 
3. The major differencing components are: 

Cll' C 13' and C 31- 

C 11 was C 3 uant ^- zed and represented by 5 information bits. 

and are each quantized and represented by 3 infor- 

mation bits. The reference frame was 2 dimensional com- 
pressed with an average of 2 bits per picture element. 
Hence, the overall data requirement is slightly over 1 bit 
per picture element. 

This method has the advantage of using only a 
partial frame memory. Here, (for the encoder) only the 
information of major components of the reference frame 
is required. Moreover, the information storage may not 
necessarily be the original data with 8-bit precision. 

This method was implemented, and it was shown to be 
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capable of reproducing a scene containing normal object 
movement without degrading the picture quality beyond 
recognition. The objective of this report is to explore 
and search for possible improvements in this technique. 
These are summarized as follows: 

1) The data compression on the reference frame 
is modified per Section 3. There, a new set of quantiza- 
tion and decoding strategy is used. 

2) This study found that it is preferrable to 
update the differencing components with new differences. 

(In contrast, in the previous LINK ABIT video study, the 
differences were always obtained between the new and the 
reference frames) . In this manner, the coding can be 
operated in true DPCM mode, and, the errors can be suc- 
cessively corrected. This is particularly desirable 

for scenes with short movement duration, i.e., containing 
objects that move and pause, it enables better response 
to stationary objects. 

3) In the previous LINKABIT video study, the d.c. 
component, C.^, is quantized to 6 bits in the reference 
frame and by 5 bits in the differencing frame. The improve 
ment using DPCM coding reduces the requirement of the d.c. 
component in the reference frame to 5 bits. Thus unless. 
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ih the differencing frames, this component can be coded by 
4 or less information bits, coding this component as 
differences does not offer special advantage. Experiments 
have shown that it is generally insufficient to use 4 in- 
formation bits to encode the frame differencing d.c. com- 
ponent. This can be explained by the fact that human eye 
is quite sensitive to object displacements, and the fact 
that the time statistical correlation between corresponding 
elements of two adjacent frames, given motion existed, is 
rather small. (In fact, most statistical correlation 
exists due to stationary objects; its advantage can only 
be extracted by motion indicators) . 

Hence, it is preferrable to transmit the d.c. component 
in the regular mode, i.e., logarithmic DPCM coding with 5 
information bits. In this manner, the requirement of frame 
memory can be further reduced. 

4.1 • Quantization of Buffer Free Frame-to-Frame Dif- 


ferencxng Technique 


Based upon the above discussion, the quantizations for 
the buffer -free frame-to-frame differencing technique is as 
follows: 

1) The reference frame is encoded as described in 
Section 3. The corresponding quantized table is given in 
Table I . 

2) The d.c. component of the differencing frames is 


quantized as in 1) . 
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3) The Hadamard component C^ 3 and C-^ are repre- 
sented by 5 information bits. Their qucmtized and trun-: 
cated value obtained in the reference frame are storage 
in the frame memory (5 bits each) . In the difference frames, 
the difference is quantized by Table II (using 3 information 
bits each) . The frame memory is updated by the quantized 
values. The total storage requirement for the frame memory 
is (for the encoder) : 

2 x (512/4) x (480/4) x 5 =153 • 6K bits 


! 1 



■ j 






The 5-bit representation of C 13 and C 31 does not have to 
cover the entire range of C^ 3 , C 3 ^ (i. e. , -255/rim/2 £ C^ 3 , 
C 3 ^'£ + 255/mn/2) . It has to cover only the absolute 
maximum change given by the 2 dimensional quantization, 
in this, according to Table I, the 5-bit representation 
linearly partitions the range: -150 to +150. 


4 . 2 Compression Experiments 

Experiments were performed on the LIM. Video Con- 
troller for k = 3, 5, 7, and 9. Correspondingly, the coding 
efficiencies are 1.0156, .9Q625, .8516, and .8187 bits per 
picture element. The results, together with the original 
A-to-D/D-to-A and the 2-dimensional compressed pictures, 
were recorded on video tape. The tape presentation con- 

' • \-L- * 

sists of processed video first run at normal speed (i.e. 

30 frames/second) then immediately followed by instant replay 
at 8 to 1 slow motion rate. The tape was recorded by a Sony 
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and Cg^: 


Quantized by 8 Levels (i.e. , 3 Bits) Each 


Outpoints Representative Value 


67 

87 

37 

52 

10 

23 

-8 


0 

-32 

- 

20 

-50 

. 

' 

41 

-75 

- 

62 


-100 


Table II. Quantization Table 
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video tape recorder, which unfortunately has very limited 
bandwidth and occasionally develop an appreciable amount 
of noise. 


From the experiments, it is evident that stationary 
objects in a scene are reproduced with reasonable quality. 

But for objects with high rate of displacements, edge distor- 
tion is clearly visible. The effect becomes highly objec- 
tionable to the observer for large value, k. The edge 
distortion is clearly caused by the lack of high order 
Hadamard components. This effect can be reduced by intro- 
ducing more differencing components. This method will in- 
crease the size of the frame memory and the average infor- 
mation bit required per picture element. Time sharing, in 
the form of updating some components at one frame and 
some others at different frames, can also improve the motion 
quality. 


Frame sharing or updating the differencing components 
was tried. The differencing strategy is given as follows: 


1st and 5th differencing frame: 

C 13' C 31 / C 12 at 3 bits each 
2nd and 6th differencing frame: 

C 13' C 31' C 21 at 3 kits eac h 
3rd and 7th differencing frame: 

C 13' C 31' C 14 at 3 kits each 
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4th and 8th differencing frame: 


C 13 f C 31' C 41 at ^ bits each. 


repeating every 9 frames. Each differencing frame uses 14 
bits. The overall average efficiency is 1 bit per pel. 

The experiment is recorded on video tape (Sample E) . 

A computer simulation program is given in Appendix B. 
From the simulations, the technique described in the above 
method provides a scheme of transmitting picture informa- 
tion with relative ease for hardware implementation. The 
degree of quality depends mainly on the factor k, the num- 

i 

ber of differencing frames, and the number of time sharing 
differencing components (the hardware complexity increases 
directly with the number of differencing components used) . 
Hardware implementation for the buffer-free frame-to~f rame 
differencing technique is summarized in the following sec- 
tion. Detailed description of individual functional blocks 
is skipped. The main points stressed are the general hard- 
ware involved and the bulk storage required for the encoder 
and the decoder. 

A possible improvement on the above technique is by 
interpolation of high order Hadamard components (those that 
were discarded in the differencing mode) between reference 
frames. This method has not been simulated. 


41 


4.3 Requirement for Hardware Implementation 

A buffer-free frame-to-frame differencing technique 
based upon the discussion of Section 4.1 can be implemented 
with relatively simple electronic hardware. The block 
diagrams of the encoder and decoder are illustrated in 
Figure 7 and 8, respectively. 

For the source encoder, due to the horizontal orienta- j 
tion of the regular TV signals, a line buffer (storage for 3 
horizontal lines + 4 samples) is required in order to per- 
form the 4x4 Hadamard transformation. This line buffer 
normally requires high read and write speed (these are 
determined by the period of the A-to-D sampling clock pulses) . 
The Hadamard transformation is performed by a fast 
serial/parallel Hadamard transformer. Since each 4x4 
Hadamard transformation requires 16 operations of add or 
subtract, sufficient bits should be reserved to avoid error 
due to truncation. 

Of all the Hadamard components, only the d.c. component 
is DPCM coded. This requires one additional data storage, data 
comparison and updating. Quantizations are carried out by 
ROM (read-only memory) table search. For components with 

i ■ 

bit- sharing, additional coding, such as bit mixing, is re- 
quired. For example, the components C 33 , C 34 , and are 
coded by 7 bits, each with 5 levels of quantization. This 
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can be done as follows: can be represented by an integer 

having value between 0 and 4. C can be represented as an 
integer of the form 5k, k = 0,1, 2, 3, or 4. C^., can be 
represented as an integer of the form 25k, k = 1, 2, 3, or 4# 
The resultant code is the sum of the integer representations 
which has range between 0 and 124 (this is readily repre- 
sented by 7 bits) . If C is the resultant coded integer, then 


C43 can be recovered by the residue of C divided by s. c 34 


Can be recovered by the residue of [C/5] (i.e., the quotient 

of C divided by 5) divided by 5. And can be recovered 
by the quotient of C divided by 25. 

The frame storage of the encoder is required to store 
the updated information of the differencing components. These 
components are updated to their nearest representative value. 
For component and C^r (each are represented by 5 bits) 

the frame memory of size 153.6 K bits is required. (If 
additional differencing components are introduced, the frame 
memory will increase correspondingly) . Speed requirement for 
the frame memory is generally very small and can be implemented 
by many low-speed, low-power devices. 

At the receiving end, the decoder requires an almost 
complete frame memory. Here, data that is transmitted as 
differences or discarded in differencing, must be preserved. 
Only the d.c. component, which is transmitted continuously, 
does not require storage. The differencing components, C, , 


and C^/ each are represented by 5 bits. The overall size 
of the frame memory is 

(512/4) x (480/4) x (32-5+2) = 445.44 K bits. 

The frame memory stores only the encoded data for those com- 
ponents discarded in the differencing frames (but not dis- 
carded in the reference frame) , and are retrieved by the 
inverse quantizer. For the differencing component, they 
are stored as the 5-bit representation of their updated 
inverse quantizer. As in the encoder, this memory can be 
implemented by many low-speed, low-power devices. 

The inverse quantized Hadamard components are fed 
the inverse Hadamard transformer. The resultant data is 
limited to the usable range of the D-to-A converter and 
stored in the 4 line buffer before it can be converted 
into video signal by the D-to-A converter and sync adder. 

The frame memory size of the encoder and decoder can 
be decreased slightly if the updated differencing components 
are .approximated by their 2- dimensional quantization (4 bits 

eaich is sufficient) . In this case, the frame memory size 

■ ! . 

for the encoder is 122.88 K bits, and that for the decoder is 
414.72 K bits. However, a quantizer is needed in the decoder 
and the appoximation using this method will not be as accu- 
rate as the previous case. 
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4.3.1 Hardware Estimations 


A general estimation of power consumption, size 
and weight for the above system, based upon the LINKABIT 
Real-Time Video Compression System, model LT8E/LT8D, is 
itemized in the following paragraphs. (Power supplies, 
their efficiencies, channel encoder, and additional data 
rate buffer are not included in the assumption) . 

4. 3.1.1 Video Source Encoder 

The video source encoder can be partitioned into 
the following submodules: 

1) Timing Generator : which generates all pertinent 

timing signals, such as sample clock pulses, control signals 
etc. 

2) Hadamard Transformer : which performs the fast 

Hadamard transformations in a "pipe- line" configuration. 

3) Line Buffe r: which stores 3 horizontal lines 

of data prior to the Hadamard transformation. 

4) Quantizer : which provides all 2-dimensional quan 

tization and quantization for the differencing components. 

5) Frame Memory : this storage in the encoder is 

used to store the information of the differencing components 
It can be implemented by high density CCD (Change Coupled 
Device) shift registers. A low power and compact version 
of size 16.384 Kbit is available for such an application. 

Ten such devices are sufficient for the encoder frame memory. 
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6) A-to-D Converter and Other Related Analog 

Circuitries : this provides proper d.c. res- 

toration and conversion of video signals into digital data 
prior to processing. 

The estimated power dissipation and number of ICs 
required are summarized in Table III. 



ESTIMATED POWER 
DISSIPATION 

ESTIMATED NUMBER 
OF IC’s USED 

Timing Generator 

7.5 watts 

I 

40 

Hadamard Transform 

18 watts 

56 

Line Buffer 

13 watts 

35 

Frame Memory 

7 watts* 

25* 

Quantizer 

10 watts 

50 

A-to-D Converter 
and Related Analog 
Circuitries 

13.5 watts 

Size: 100 cubic 
inches 

Weight: 2.6 lbs** 


* This includes data control logic. 

**This is based upon the A-to-D converter manufactured by 
Mirco Consultants, model AN-DI-802 RAD-B, which has a 
size of 3.7" x 2" x 9.5” and weights 2.2 lbs. 


TABLE III 
















Using the assumption of packaging density of 
1,5 IC ' s/cubic inch 
and 30 IC s/pound, 

the estimated power dissipation, weight and size are: 
Estimated Power Dissipation: 70 watts 

Estimated Weight: 10 lbs 

Estimated Size: 240 cubic inches 

4. 3.1. 2 Video Source Decoder 

The Video Source Decoder can be partitioned into 
the following: 

1) Timing Generator, 

2) Hadamard Transform, 

3) 4-Line Buffer, 

4) Inverse Quantizer, 

5 ) Frame Memory , 

6) D-to-A Converter and Sync Adder. 

The estimated power dissipation and number of IC's required 
are summarized in Table IV. 

The estimated power dissipation, weight, and size are: 
Estimated Power Dissipation: 83 watts 

' - i | 

Estimated Size: 190 cubic inches 

Estimated Weight: 8.5 lbs. 
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ESTIMATED POWER 
DISSIPATION 


ESTIMATED NUMBER 
OF IC'S 


Timing Generator 


7.5 watts 


30 


Hadamard Transform 


18 watts 


56 


4-Line Buffer 


15 watts 


36 


Frame Memory 


20 watts 


50 


Inverse Quantizer 


15 watts 


50 


D-to-A Converter 
& Sync Adder 


7 . 5 watts 


Size: 40 cubic ins. 
Weight: 1 lb. 

(Self-enclosed in 
box) 
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A P P E N-D IX A 

The objective of this Appendix is to sketch briefly 

the proof of the inequality 
m n 

£ £ 

i=l j = l 


| c i j J - I C 1:L I + I C ;L1 *-255^mn/2 1 < 255/inn/2 


It should be helpful to observe the case when m = 1 and 
n = 2. In this case, the.Hadamard basis is 



— ( 1 , 1 ) 

/2 



— ( 1 , - 1 ) • 
/2 


The original domain of the 1x2 subpicture and the trans- 

4 

formed coordinates are illustrated in Figure Al. Since the 

_y/ 

transformed components must lie within the square of the 
original signal domain, then the components, according to 
the figure, can be easily verified to satisfy the follow- 
ing .inequality 

|C, , - 255/ /2 | - + }C, J < 255//2 (kth nearest 
11 integer) 


which is the above inequality. 

k 

For the arbitrary case mn = 2 , k >_ 1, it can be 
readily proved by first coordinate translation of the 
-original subpicture, Y 


x 


11 


k ln 


Y 


ml 


x 


mn 
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into 


x' . . . x' 

■■ (: .) 

' X n . . . X ' 

ml mn 

when x ' . = x . . - 255^/2. 

lj 

x' . has range -255v / 'mn/2 to +255v / ^rm/2. (Note: all 

ID 

components, evcept the d.c. component, of the Hadamard trans 
form of are identical to those of Y) . Then apply the 
fact that the inverse Hadamard transform in the translated 
coordinates has peak value (i.e., max ! x ^j|) equal to the 
sum of the absolute values of the individual Hadamard com- 
ponents (which can be shown by induction on k) . 

Remark: Source coding without using the inequality 

(3. 14) may result in inverse" transformed signals lying 
outside the domain of the original signal. Using the in- 
equality helps to confine the occurrence of video signals 
of the ensemble of recognizable picture to more likely area 
of occurrence. 
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APPENDIX B 


F? 

W 


t 


The computer simulation for the buffer- free frame- to- 
frame differencing technique is simulated by the LIM Video 
Controller. The Fortran program is listed in the following 
pages. Most of the subroutines used are explained in the 
LIM Video Controller Operations Manual. The following is 
a list of subroutines used in the program but were not in- 
cluded in the LIM Video Controller Operations Manual. 

1) MEAN: Calculate the arithmetic mean of a set of 

numbers. (Note: in the simulation, the average of the d.c. 

Component over a horizontal line group is used. This is used 
because of hardware implementation of the LIM Videc Controller 
where substantial sampling occurs in the blanking and sync 
region. 

2) DBNQT , DLNQT : Calculate, successively a set of num- 

bers, their successive differences quantized and updated by 
inverse quantizations. This simulates the DPCM coding. 

DBNQT uses logarithmic search technique, while DLNQT uses 
linear search technique. The quantization outpoints and in- 
verse representative values are specified as pointers in the 
arguments of these subroutines. 

3) DIFF and RENEW: These subroutines substract and add, 

respectively, two sets of numbers. They are used to calculate 
differences and perform updating for the differencing components. 
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4)1 BREAD, ZWRIT, ZTEST and ZINIT: are special fast 

digital disk routines. 

Titles are prewritten on the disk with name "TITLE." 

The program allows simulation of time-sharing dif- 
ferencing components. The differencing information is stored 
in the array, NZX, which indicates which components are to 
be differenced and provides pointers for the quantization and 
inverse quantization allocated for the component. 

In the video simulations, a grey vertical band is 
visible to the right side of the picture. This is the glitch 
filter override, which overrides the switching spikes caused 
by the video disk recorder. 
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/ FOR 


♦LIST ALL 

C S-CD-LS-IQL-. 

*I0CS(2501 READER) 

*IOCS( KEYBOARD, TYPEWRITER) 


STNO.C 


FORTRAN SOURCE STATEMENTS 


I DEMTFCN **COMPILER MES 



INTEGER J D 1 (4162) , J D 2 (41 6 2) ,B F 1 (204 9,2 ) ,B F2 ( 2049,2 ) 
INTEGER 0 L2 (3?) , IQ L2 (3? ) ,Q1 (16), 101 (16),Q3(9),IQ3(9) 

UiJ E G ER._Q.4 (7),IQ4„m,Q 5 CAXrJLQl 111 " 

INTEGER NZX<15,4),DQ(16,3),IDG(16,3) 

INTEGER CHAR(16,16),DIGIT(3),FNU!'U2),EOF,BOF 

EQUIVALENCE — C J D 1 (.1 ) ,B F I ( 1,A)),AJD2X1) ,B F ?.(1,1)> 

DEFINE FILE 1 (730 ,320 , U ,K ) 

DEFINE FILE 2 < 1 30 ,52 0 , U , K ) 

DEF INE—.EJL LE_3 ( 4 0, 25 6, IF, K) 

DEFINE FILE 4(130,320,0,0 

DATA FNUM/1 ,4/,JDl (41 62) , J D2 (4162) /2*4160/ 

D AT A — D IG IT 744 0,420, 4007, L EV-tTSQT^L E V 2/100/ - 

CALL ZINIT(FNUN,2,IFG) 

READ (2*1) B F 1 

CAL 1 ZW R-I-TA 2 ,2 ,-J D_1X4 1.62JJ 

READ (8,100) E0F,NLG,NFRM,IFG,M1 
WRITE (1,200) E0F,NLG,NFRM,IFG,M1 


396 Ml =0 


GOTO 397 
95 _W 1=2 


397 n0F=E0F-4*NLG+4 
CALL HMODE (M3,M4) 

D.0_5A6 _L= 1,-4 1 

READ (8,100) (NZX(J,I),J=W15) 
IF (IFG) 566,566,588 


iHMmuvAimi 


OF POOE QUALITY 







310.X FLjOR T_ 


588 WRITE (1/200 (NZ X ( J / I) =1 / 1 5 

S66 COMXJ-MUE 1 , ' 

READ (8/100) 0.L2 
READ ( 8/ 1 00) 1 Q L2 

— R E A D — (S, ID Q X_Q 1 __ 

READ (8/100) 101 
READ (8/100) 0.3 

; — READ — (£ / 1 00) — I Q3 

CALL ZTEST (2/IX/1 ) 

READ (2' 14) 8 f 1 

— - ■ C ALL- 7 U R I T C2 /-1./-J.D-1-C4-1420 ) 

READ (8/100) R4 
READ (8/100) IQ4 

READ -C 8/ .1 <10)_ U 5 

READ (8/100) IQ 5 
IF (IFG) 411/411/413 

■44-3- W RITE— (.1 /333) 

WRITE (1/200) QL2 
WRITE (1/250) IQL2 

MS I T E — ( 1. / 333) 1 

WRITE (1/200) Ql 
WRITE (1/250) I Q 1 

_.W R I T EX 1 / 3 3 3 ) , 

WRITE (1/200) Q 3 
WRITE (1/250) IQ3 

WRITE-C1 /333> 

WRITE (1/200) Q4 . 

WRITE (1 /250) IQ 4 

WRITE— (1/33 3) : 

WRITE. (1/200) Q 5 
WRITE (1/250) IQ5 

-4J-1— W R-I T E— - (1/353 ) 

DO 197 1=1/3 

READ (8/100) (DG(J/I)/J=1/16> 

READ — (8 /.I 00.) — (.1 DQ (J *X) / J-= L/-16-X- 

IF (IFG) 197/197/194 
194 WRITE (1/200) (DQ(J/I)/J=1 /16) 


— E_J4—X_S X DENTf CM . *» C Ottt-MLER- : 


Read and Output to the Console. The Quantization 
Outpoints and their Inverse Representative Values. 


ORIGINAL ;; 

OF POOL QLALITA 




F CL- 


WRITE (1*250) (IDQ(J*1)*J-1*16) j 

WRLT-E — C1.*_33 3 ) T — 

197 CONTINUE 

CALL ZTEST^IFG*!)^ 

DO— BO S-I E.R>1=L* N ERE4 !! 

IDFX=MOD(IFRM-1*9) 

I LN-EOF 

CA L L_ F-A LOD_(_l./.I LN.*.Q.) 

I R F =2 

DO 800 1LG=1 *NLG 

CALI ZJ£ S T_(_l.rXF.G /_!_) 1 

IF (IDFX) 355*355*356 
356 CALL ZREADd *ILG*JD2(4162>) 

C ALI Z-TE S U(4-*_L£ G * -U 

355 DO 269 I P=1 *2 

RF1 (2049*IP>=2048 \ , 

CALI R D 4 LN-tU# BF4.(-2Q49*-I.&)4 

IF (ILG-NLG) 251*248*248 
248 IF (IP-1) 251*251*259 

-254— CA L i FJLT-QD X 1R£*.I LN/4) — 1 

IF CIRF-1) 258*258/256 
258 IRF=2 

G 01.0-25.9 

256 I R F = 1 

ILM=ILN-4 

-25-9 C A L I S 3 H ft.( S.FJL (20 4.3* IP)-* a F_lX1.536-r L® ).* 5 1.2*420)— 

CALL SHM2(BF1 (2048* IP) *RF 1(1 024 * 3 P ) * 1 024 ) 

CALL WEAN (125*BFK2048*IP)*BF1(1*IP)*0) 

CA LL- DBM QXC.3 L 2 ( 3 2) * 3 2. * I Q-L2-C 3 2) *-_3 F_1 .( 204 8*.X EX*. 

IF (IDFX) 347*347*337 
347 CONTINUE 

C A LL-LNQ T-ta.4- C.7.) * 7.* I Q 4 ( 7 ) *.B F 1.C2Q 47-,.LP_)*_1 2-8X 

CALL LNGT (Gl (16)*15*IQl ( 1 6 ) *B F 1 ( 2046* 1 P ) * 1 28) 
CALL LNQT (03 (9) *9*103 (9) *0F1 (2045 *1 P) * 128 ) 

C All l ^ QT(Q4 ( 7-)-/X/-LO 4 CZX* D F4X153A* I g-).XL2.8) 

CALL SETVLC BF 1 (1 535* IP) * 1 28*0) 

CALL SETVLC B F 1 (1 534* I P ) * 1 28*0 ) 


Starts Processing 

K + 1 (To change the 
is replaced by IDFX = 


value of K 
= mod' [ifrm 


this card 
HIV k + 1] ) 


Dmponenl 


OP POOP QLiALl. 







E — 0— R_-T__R_ft_ -N S__ a_UL_R__C_E S . .T S_ X_ JE. J4_F— M— T— S 

CALL SET VL ( BF 1 ( 1 533/ IP ) xl 28 x0 ) 

-CAL I — LH.Q T..CQ.1. LI 6 X* 1 5.x IQ 1 C1.6.).x B FLtL024xLPlxX2.81 i 

CALL SETVLC RF 1 <1 023x TP) x 1 2BxP) I Quan 

CALL LNQTCQ5(5)x5xIQ5(5)xBF1(10?2xIP)x128) f Refe 

CALi L.NQT..CQ5 ( 5.) x 5.x 1G5 (.5 )xBF.1 (1U2.1_x.I Bl x.1 2iD __ 

CALL LNUT(Q3(9)x9xIQ3C9)xRFl ( 5 1 2 x I P ) x 1 28) . 

CALL SETVLC BF1 (51 1 x IP) xl 28xU) 

CA L I AN QT.CQ.5.C Six 5 x. LCL5-C5.lxS.Ll 1510 x I S)-xX2.SX 

CALL SETVLC B FI C509xl P) xl 28x0 ) J 

GOTO 269 

COMT-I HUE ll . 

CALL XFER (BF1 C2048xIP)x3F2(2048xIP)x128) ^ 

I Y =r*!0 0 (I D F Xx 4) +1 

D. 0.„3 0 5— I.=1 x.1.5 . : 

IX=NZX CIxTY) 

I FG=I / 4 

IXG-X 4-1 FG) *5 12 1 

I FG = I. FG-flOD (1x4) 

IF (IX) 80 5 *,>805 x 8 86 

88 6 CALL-D IE F-C8 E.I Cl F Sx IPl-xaF-2..aj-Gx lPlx-1281 

’CALL LNQT ( DQ ( 1 6x I X ) xD Q ( 1 x IX ) x ID© ( 1 6 xl X ) xB FI (IFGxI P)x1 28) 
C ALL R FNE W (BF2 ( I FG x IP) xB FI ( IFGxIP) xl 28) j 

^—CONTINUE 

269 CONTINUE 

CALL ZTEST (1 xIFGxI ) 

LF Cl D F X 1_ J?-9_Lx 99-lx-992 I 1 

991 CALL ZWRIT (1 xILGx JD1 (4162)) > Store in Digital Disk 

GOTO 800 j 

992 CA L I Z W R 1X1X# 1 LG x J D 2 C4XS21 ijL 

800 CONTINUE 

CALL RHODE (HlxM2) 

LLN =E0T - - c 


Quantization and Inverse Lookup for 
Reference Frames 


Quantizing and Inverse Lookup 
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CALL BLANK (1xE0Fx2xB0F) 
CALL 7 TEST (1 xIFGxI) 

-CAL t ZR.EAJ)_(J-x.1-xJL DXJL4 162- 

CALL BLANK(2xE0Fx1 x80F) 


OF POOR ODAl 



IDK=2 

D 0-7 9 7—M T=X,N LG- - : 

CALL ZTEST (1,1 FG,1 ) 

IF (TLG-NLG) 166,166,985 
446—11 — (I D K - 1 >—98 4,986 ,.9 86 • 

984 CALL ZREAD (1 ,ILG,JDl (41 62) ) . 

GOTO 935 

-9-&6_ C ALI Z REAQ11, 1 LG,_J D2-C4.1 6-2J J_ 1 . 

985 DO 757 I F ’1 = 1 ,2 

IF (IDK-1) 381,381,382 

382 SF.1 (2049, 1 HO =2 048. 

CALL SBHM (8 FI (2048, I F PI) , B F 1 ( 1 5 36, I FM> ,512,128) 
CALL SHH2 (BF1 (2048, I FM) ,8 FI <1 024 , I FH) , 1 024) 

-JCA LI — LMRN4 CBF1 (2048, I FM ) , 2.0 4 8,-0^ 409 a) 

CALL WR4LN(G,GF1 (2049, IFM>) 

GOTO 757 

-381— B F 2- ( 2 0 49 , I F .» ) = 2 0 48 ; „ 

CALL SBHM(BF2 (2048, IFM),BF2(1 536, IFM), 512, 128) 
CALL SHM2 (8F2 (2048,1 FH) ,B F2 (1 024 , I FH). ,1 024 ) 

,C A L L-LN R N 4 ( B F 2 ( 2 0 4 3,1 -F-M )-,-2 04 8 , 0,4090-) 

CALL W R 4 L N (0,BF2(2049,IFM)) 

757 CALL FDTOA <IFH,ILN,1) 

1LG-I LG + 1 — 

IDK = F!OD (1DK,2)+1 
797 !LN=ILN-4 

— CALL - ZTEST (-1 , I-FG,-1-) - 

DO 681 I X =1 ,2 
DO 681 IGP = 1 ,2048 

48-1— O F 1 (-1 GP, I X ) ~LE V2 : 

C ALL UR4LM (0,BF1 (2049,1 ) ) 

1P=ILN 

-MT - I L-M - 1 2 — : 

DO 466 IFM=1,2 


GOTO 


(771 ,772), I FH 
(T, I-F-GVD- 


CALL 

GOTO 


ZREAD (2,1 , J D 1 (4162)) 
773 



"V 


. Inverse Hadamard Transform, Limit to 
> Usable Range, and Record on Video Disk 
Record 


J 




- O RIGINAL PACtF, IS 


OF POOR QUALITY 



:» 


r< 


772 CALL CTEST C2,IFG,1 ) 

CALI Z-R£ AJU2.*2L* J.D.2 ( 416?>; 

773 CONTINUE 

CALL FDTOACIFM,MT-3,0) ■ 

—CAL I FA CO A-Cl-FALrilT^^Ili 

CALL FDTOA (I FH*HT*0> 

466 CALL FDTOA (IFM^ILN^I) 


; 

CALL /_TEST<2*1FG*1 ) 
DO 682 I LG =1 
IY-MOD l T„10>+1 


1 = 1/10 

R t AD (3 * I Y) CHAR 


DO 364 JL.\I = 1,4 
DO S64 I FM=1 *2 

IX=51 2 * C 5 -1 LN ) -D LGI T-CILG) 

DO 865 MT=1,16 
I F 6 = C HAR (MT^IRX) 

CALL—SGN CIFjGJ 

BF2 (I X/IF?0=IFG*LEV2+LEV1 
I FG-CHAR ( MT^IRX+S) 

CALL...SGN C LEG) 

8 FI (I X*IFM)=I FG*LEV2+LEV1 
865 IX=IX-1 

864 m y.-i Rxt 1 

682 CONTINUE 
I LN = I P-4 

ID Q - 69.1. I Fa!~l/ ? , 

GF1 (2049r I F M) =2048 
B F2 (2049* I FM) =2048 

CAL! — UR4 Ui CU*BTlX20AA*XEJtt->_ 

CALL FDTOA(IFM*ILN*0) 

CALL WR4LN C0*BF2 C2049*IFM) ) 

LI — F-D-TTLACI-LM ,1 lN-4^0) 

MT =M1 
Ml =1X2 




» 


I 


M2-MT 

3 

M3=M4 
H4=MT 

3±_1_ 

CALL MOVE (9#1 *MT) I 

CALL MOV E (1 r MT ) j 

&D8 CONTINUE J 

CALL MOVE (0*0-2) 

WRITE (1*444) 

CALI EXIT J 

100 FORMAT (815) 

200 FORMAT (5X,8I10> 

H- F.O B S A.T.C8 LAID ■- 

333 FORMAT (»****») 

444 Format ( * * end processing**) 


Go to Next 


ARIA9LE ALLOCATIONS 


J 01 ( 1*2 

LQL2.(J_*2_ 

G 4 ( I* 2 
D Gt ( I * 2 

£Of,LL* 2 !. 

NFRMC 1*2 
I (1*2 


I DK ( I* 2 
I R X ( I * 2 


>=105D-001C 
_Jl-2QE.1-2C.C2- 
) = 2 1 1 A- 2 1 1 4 
>=2197-2168 

)_= 22 CD. 

>=2202 

>=2207 

D.C 

>=2261 

>=2266 


B FI ( 1*2 > 

Q.1 Cl *2 >. 

IQ4 ( I* 2 > 

109(1*2 > 

BO F.(I*2 ) 

MIC 1*2 > 

J ( I *2 > 

— LR FX1 * 2 >. 

HT ( 1*2 > 

LEVI (1*2 ) 


=1 05D-005C 
* 2011-20 E2- 
= 2121-2110 
=21 C7-2198 

=. 22 .CE 

=2203 

=2203 

=22QD 

=22E2 
= 22E7 


STATEMENT ALLOCATIONS 


_4 .1.3 =24 3.6. 
259=2580 
984 =28 E A 


200=2306 
_4J 1=2.49. CL 
347=2618 
986=23F8 


250=230 A 
_L94=2„4CC_ 
337 =27 9D 
985=2904 


333=2300 
J 97-= 25 J 2 
886=27 EF 
382=290E. 
_8.6i4-=2B_l.D_ 


IFEATURES SUPPORTED 



444=2311 

_35l6=253A_ 

8Q5=?36A 

381=2972 


396=236A 
_3 5 5 =25 4-EL 
269=2873 
757=29D4 


395=2 37C 
-248=256C_ 
991=2885 
797=29F0 


397=23«C 
251=257? 
992=2893 
681 =2 A 0 C 
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